Stretching semiflexible filaments with quenched disorder 
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We study the effect of quenched randomness in the arc-length dependent spontaneous curvature 
of a wormlike chain under tension. In the weakly bending approximation in two dimensions, we 
obtain analytic results for the force-elongation curve and the width of transverse fluctuations. We 
compare quenched and annealed disorder and conclude that the former cannot always be reduced to 
a simple change in the stiffness of the pure system. We also discuss the effect of a random transverse 
force on the stretching response of a wormlike chain without spontaneous curvature. 
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Introduction. — The wormlike chain (WLC) model 
of semiflexible polymers treats the macromolecules as 
one-dimensional locally inextensible curves with bending 
rigidity [iHl]- Despite its simplicity and its initial suc- 
cess with the experiments on stretching of ds-DNA 
the inherent complexity in the microscopic architecture 
of biological filaments invites the search for more realistic 
models. One key feature of many biopolymers, such as 
DNA or denatured proteins, which goes beyond the clas- 
sical WLC is their heterogeneity. Their local architecture 
(e.g., spontaneous curvature) 7] or their local elasticity 
(e.g., their bending rigidity) vary along the polymer con- 
tour. At large scales, this heterogeneity may statistically 
behave as a random variable. It may be fixed in time 
such as that related to the different base pair sequences 
in DNA. This randomness is treated as quenched dis- 
order. If the random inhomogeneity undergoes thermal 
fluctuations, it is treated as annealed disorder. An ex- 
ample of the latter is the reversible binding of proteins 
to DNA capable of inducing spontaneous (intrinsic) cur- 
vature which has been studied in [8|. 

The first studies of the effect of random base pair se- 
quences on the macroscopic statistical properties of DNA 
focused on the related random bends (spontaneous cur- 
vature) and treated them on equal footing as the ther- 
mally excited bends, thus considering the annealed disor- 
der case . Bcnsimon et al. considered the effect of 
random angles in the Kartky-Porod chain and using nu- 
merical transfer matrix and Monte-Carlo methods found 
that, as far as stretching is concerned, the disordered 
chain behaves as a homogeneous chain with a renormal- 
ized persistence length [12j. P. Nelson calculated the ef- 
fect of weak disorder on the entropic elasticity of a flexi- 
ble rod with random kinks and found that the sole effect 
was a reduction in the apparent bending stiffness with 
the twist stiffness remaining unchanged [13| . A different 
way of modeling disorder has been followed by Debnath 
and Cherayil [14j who considered the chain under ten- 
sion as consisting of random A-B blocks with different 
elastic constants. The blocks were treated as "Gaussian" 
semiflexible polymers (without the local inextensibility 
constraint). Muhuri and Rao [IH studied finite size ef- 



fects in the Kratky-Porod chain with a random sequence 
of stiffness constants. 

In a recent paper [lj|, we studied the response of a 
weakly bending WLC with arc-length dependent sponta- 
neous curvature to a stretching force applied at its ends. 
We specifically considered the case of sinusoidally varying 
spontaneous curvature which allows us to treat the gen- 
eral case via Fourier transformation. This simple and an- 
alytically tractable model appears to be particularly well 
suited for the study of quenched disorder in the spon- 
taneous curvature of a filament under tension. This is 
the subject of the present paper. A similar calculation 
with a different method (which is in principle valid) was 
attempted in Ref. [17], but a mistake in the way the 
thermodynamic limit was taken led to the incorrect con- 
clusion that sequence disorder has no effect on elasticity. 
Here, we use the replica trick to calculate the effect of 
uncorrelated quenched disorder and we obtain the force- 
extension relationsip as well as force-transverse- width re- 
lationship. We compare our results with those obtained 
in the context of annealed disorder. Although, as we 
show, the two cases (quenched and annealed) are indis- 
tinguishable in the limit of weak disorder, our results for 
quenched disorder hold for arbitrary strength. 

We also study the adaptation to the WLC of the clas- 
sical random-force model which Larkin had introduced 
in the context of the Abrikosov lattice 18]. As with 
quenched disorder in the spontaneous curvature in the 
weakly bending approximation, the coupling of a random 
force is linear and this makes its analytical treatment 
straightforward. The effect of quenched random forces 
on the WLC conformations has recently been analyzed 
in [ijj]. Here, we deal with their effect on the stretch- 
ing response in parallel with our study of the random 
spontaneous curvature. 

Model. — We consider a WLC fluctuating in two di- 
mensions under a strong stretching force applied at its 
end-points. The restriction for the chain to remain con- 
fined to a single plane is made to allow concise analytical 
treatment. It certainly applies to stretching experiments 
on a surface [13] • Although it is known the difference 
between 2d and 3d, when the curvature of a freely bend- 
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ing WLC is concerned, to be non-trivial 2l|, here we 
stay within the weakly bending approximation where it 
can be shown, using the Monge gauge [22}, that the two 
transverse directions of a (1 + 2)-dimensional chain de- 
couple. Our recovery, for a special case of disorder, of 
the main result of Ref. [l3} which was based on a 3d 
model further attests to the validity of the main features 
of our results beyond 2d. Quenched disorder along the 
polymer contour is represented by a random distribu- 
tion of the arc-length dependent spontaneous curvature. 
The strong stretching force allows us to use the weakly 
bending approximation which significantly simplifies the 
analytical treatment of the filament response as shown in 
Ref. [lfij]. The WLC conformations are parametrized by 
the displacement y(s), transverse to the direction of the 
pulling force (x), as a function of the arc-length position 
s. The elastic energy functional is given by 
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where k is the bending rigidity related to the persistence 
length via n = ^L p k B T, f is the applied stretching force, 
L is the total contour length, and c(s) is the signed spon- 
taneous curvature. The probability distribution of c(s) 
is given by 



"P[c(s)] ~ exp 



(2) 



and is completely characterized by c(s) = and 
c(s)c(s') — A c S(s — s'), where the overbar denotes an 
average over the disorder. The delta function reflects the 
assumption of uncorrelated disorder along the polymer 
backbone. 

Experimentally observable quantities of long polymers 
with quenched disorder are expected to be self-averaging 
and are calculated by averaging the free energy over 
the distribution of disorder [23]. This requires averag- 
ing the logarithm of the partition function, a goal which 
is achieved by employing the standard replica trick [24| . 
Averaging the partition function itself over disorder cor- 
responds to the annealed system which has been studied 
111 Ref. [HI where it is called the "disorder-first system." 

Besides disorder which is related to the spontaneous 
curvature, in this paper, we also consider the effect of 
a random transverse force on the force-extension rela- 
tionship of a WLC without spontaneous curvature. The 
clastic energy functional of such a system is given by 



H g [y(s)} = 



ds 



( d 2 y(s) y 
\ ds 2 J 
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dsg{s)y(s) 



f*(T)'-" 



(3) 



As usual, the distribution of the random force g(s) is as- 
sumed to be completely determined by g(s) = and 
g{s)g(s') = A g S(s — s'). The random coupling term 
in Eq. can be viewed as expressing the interac- 

tion of a stretched semiflexible polyampholyte with a 
quenched linear charge density g(s)/E and a uniform 
transverse field E. An alternative interpretation of Eq. 
(J3j> views g(s) as representing a random interaction with 
the crowded environment in which the filament fluctu- 
ates. Strictly speaking, quenched disorder in the em- 
bedding medium (environment) should entail a term 

Jo dsg(x{s))y(s), where x(s) = s - \ f° ds' (d s ,y{s')) 2 , 
instead of the one written above, but g(s) is expected to 
capture the leading-order behavior in the weakly bending 
limit. 

Replica Trick. — The standard method to deal with 
the quenched disorder average involves averaging over n 
identical non-interacting copies (replicas) of the system 
24]. In the end, the replica limit n — > is taken. The 
free energy is calculated via hiZ = lim„^o(^" ~ l)/ n i 
where Z is the partition function of the system with a 
certain realization of the disorder. Z n can be expressed 
as Z^ = exp(--H (rep) /k B T). 

For the random-spontaneous-curvature system de- 
scribed by Eq. ([1]), the replica "Hamiltonian" takes the 
form 



K 



(rep) 



«=i 

A r K 2 



2k B T(l + nKA c /k B T) 



a,b—l 



where we have omitted constant terms and the subscripts 
a, b label the replicas of the original system. In Ref. 
[16j , we have shown that in the strong stretching regime, 
where / 3> k/L 2 , the details of the boundary conditions 
become irrelevant. Therefore we can assume hinged- 
hinged boundary conditions and a chain with vanishing 
curvature at its end-points. In this case, we can use the 
Fourier decomposition 



y a (s) = £ sin((?mS ) 



(5) 



where q m = tttti/L is the wavenumber of the correspond- 
ing m-mode, and express the replica "Hamiltionian" of 
Eq. (j4|) in a quadratic matrix form. Correlators can 
be calculated using the Sherman-Morrison formula from 
linear algebra. In the replica limit (n — > 0), we get 



(4 m) 4 m) ) = 



2k B T 



2A r n 2 



L 9m( K <?m + /) L ( K <l?n + I)' 



,(6) 



where 1 is an n x n matrix with all of its elements equal 
to 1. Using the matrix determinant lemma and taking 
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the replica limit we also calculate the disorder averaged 
free energy and obtain 



G c 



k B T 

2 ' •• /.'/; 

m— 1 

/c B T («^ + /) 



(7) 



For the random- force system described by Eq. ([3]), 
the corresponding replica "Hamiltonian" is given by the 
analoguous procedure: 



01 (rep) 

ILg 



l 

2./o 



ds 



n 



a=l 
i-L 



d 2 y a (s) 
ds 2 



dy a {s) 



2k B T 



/ ds V y a {s)y b {s) 

J ° a,b=l 



(8) 



Using similar calculations as in the case of random spon- 
taneous curvature, we obtain the correlators 



<4 m) 4 m) > = 



2k B T 



2A fi 



^9m(«9m + /) L ltn( K Q?n + fY 

and the disorder averaged free energy 
k B T 



lab ,(9) 



Gg — 



m— 1 

1 A, 



k B T q 2 m (Kq 2 m + /) 



(10) 



Force-Extension Relationship. — The average pro- 
jected length of the filament in the direction of the 
stretching force is given by 



(x{L))=L-\ 



di 



( dy{s) 
\ ds 



(ii) 

(12) 
(13) 

In the strong stretching regime, defined by / 3> f cr ^ /l, 
where f cr = k/L 2 and /l = k/L 2 (the shorthand nota- 
tion defining the characteristic force scales) , we obtain 



where 



and 



i n 

((AM)2) = Hm-^(4-)4" 
n— >o n L — J 



{xiL)) -l = -l(^) 1/2 -U,L(lf) 



fjA 1 / 2 

L 2 V / ) 8 c V / 



(14) 



Equation (|14|) is the central result of our paper. The 
first term in the right-hand side (rhs) is the well-known 
expression associated with the ironing out of the thermal 



undulations [25|. The second term comes from straight- 
ening the quenched random undulations related to the 
spontaneous curvature. This result implies that, as far 
as the force-extension response is concerned, uncorre- 
cted quenched disorder in the spontaneous curvature 
acts as an effective linear increase in the temperature 
given by k B ST = kA c /2 (using L p — 2n/k B T). We also 
point out that for given bending rigidity and disorder 
strength, the effect of disorder on the force-extension re- 
lationship is independent of the polymer size L. This 
result holds for arbitrary strength of disorder, provided 
that the strong stretching and weakly bending assump- 
tions are fulfilled. Our result does not contradict the 



claim by Marko and Siggia 25] that disorder related to 
intrinsic bends with large radius of curvature compared 
to the persistence length do not alter the large / limit. 
In fact, any arc-length dependent spontanous curvature 
whose Fourier spectrum has bounded amplitude does not 
alter this limit. As we have shown in [16j, the response 
of a weakly bending WLC with arbitrary (but bounded) 
spontaneous curvature is the superposition of the re- 
sponses corresponding to sinusoidally varying modes of 
curvature. The latter contribute terms which scale as 
~ f~ 2 and become negligible compared with the ther- 
mal term ~ f^ 1 ^ 2 as the stretching force increases. The 
crucial difference with the case of uncorrelated disorder 
is that the latter contains bends of unbounded sharp- 
ness in a similar fashion as the thermal excitations do. 
That is the physical reason why the disorder term in Eq. 
(TT41 has the same form as the thermal one. Of course, 
uncorrelated disorder with zero correlation length is a 
mathematical abstraction. In real systems, there will 
be a maximum wavenumber, <7 max , in the spectrum of 
random undulations of the spontaneous curvature. For 
/ ^ K( ?max 01' disorder approximately behaves as un- 
correlated and our Eq. (IT41 is expected to hold. For 
/ ^> ft(?m ax , the Marko-Siggia (thermal) limit for large / 
will prevail. 

An alternative way to obtain the force-extension re- 
lationship is by taking the derivative of the disorder- 
averaged free energy given in Eq. ([7]) with respect to the 
stretching force: dfG c = {{x{L)) — L)/L . As expected, 
both ways yield exactly the same result. 

In Ref. [17], it is shown that the effect of annealed 
disorder in the spontaneous curvature amounts to the 
replacement of the original bending rigidity k by 

K e ff — K /(l + A c K/k B T) . This implies that for weak 
disorder, defined by A c <C k B T/n, the approximate lin- 
ear relation holds K e ff ~ k(1 — A c n/k B T) and we im- 
mediately recover the response described by Eq. (fT4"]) . 
In this regime, quenched disorder has the same effect as 
annealed disorder. This is the regime treated in Ref. 
[l3j ]. The difference between the two types of disorder 
becomes significant in the case of strong disorder, where 
K e ff ss k B T/A c . In the quenched case, the second term 
in the rhs of Eq. (Tl4]) dominates the response. In the 
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annealed case, we would have gotten — jA c L{Jl/ f) 1 / 2 
instead which differs by a factor of 2. A more signif- 
icant qualitative difference between quenched and an- 
nealed disorder can be deduced from the form of Eq. 
©. Except for the weak regime, the effect of quenched 
disorder cannot be reduced to a simple renormalization 
of the bending rigidity. 

The force-extension relationship for the random-force 
system of Eq. ([3]) can be calculated in a similar fash- 
ion, taking into account Eqs. ((9|) or (flTjj) . For strong 
stretching, / >• f cr > Jl, we obtain 

Wf) , l(fcr\^ lA g L 

We notice that as the stretching force increases, it quickly 
irons out the bends caused by the random force and the 
response is dominated by the classical thermal undula- 
tions. On the other hand, for given stretching force, 
bending rigidity, temperature, and disorder strength, the 
effect of disorder grows linearly with the contour length 
L and becomes dominant in the thermodynamic limit (/ 
needs to increase accordingly in order to stay within the 
weakly bending approximation). This is an interesting 
similarity with the random-force-induced destruction of 
long-range order in the Abrikosov lattice according to the 
Larkin model [181 . C omparing Eq. ([151) with the result 
obtained in Ref. [16| for the force-extension relationship 
of a WLC with sinusoidally varying spontaneous curva- 
ture along the polymer contour, we see that the random 
force effectively acts as spontaneous curvature of that 
type with amplitude c e // = {AgL/Zn 2 ) 1 ^ . 

Transverse Fluctuations. — The shape of transverse 
fluctuations of a stretched filament, ((y(s)) 2 ), is a useful 
diagnostic tool of its elasticity distinct from its extension 
in the direction of the pulling force, (x(L)) [26| . Using 
the correlators of Eq. ([5]) , we calculate the width of the 
transverse fluctuations at the mid-point (s = L/2), for 
/ 3> /l, and obtain 

JT TT27- LkB Tl A c k*/ 2 1 

For strong disorder, the second term in the rhs of the 
above equation, which scales as ~ /~ 3 / 2 , can dominate 
over a range of forces before it is overtaken by the ther- 
mal term which scales as ~ / ■ In the annealed case 
the disorder-related term is absent, whereas in the case 
of weak disorder, A c <C ksT/n, it is negligible. We 
point out that strong quenched disorder in the sponta- 
neous curvature causes a qualitatively different behavior 
in the response of transverse fluctuations which cannot 
be simply reduced to a renormalization of the persis- 
tence length of the undisordered wormlike chain. We 
see that the claim by Bensimon et al. [13] that "a ran- 
dom Kratky-Porod chain with a certain type of disorder 
is well approximated by a pure chain with an effective 



elastic constant" is not totally correct. It only applies to 
the force-elongation curve and cannot be extended to the 
transverse fluctuations of a strongly disordered WLC. 

For the random-force system, the width of the trans- 
verse fluctuations at the mid-point (s = L/2), for / 3> 
f L , is given by 

77~ L x \2\ Lk B T 1 A a L 3 1 , s 

((*-=2>) ) = ^—f + ^p- 

The random-force contribution scales as a higher power 
of the filament length L and thus is relevant in the the 
long chain limit. 

Conclusions. — We studied how the uncorrelated dis- 
order in the spontaneous curvature of a weakly bend- 
ing WLC affects its response under strong tension. The 
force-elongation curve for quenched disorder is identical 
to that of a pure system at a higher temperature. The 
quenched case is identical to the annealed case in the 
limit of weak disorder where its effect can be reduced to a 
decrease in the bending stiffness. This is no longer true as 
the disorder becomes stronger and our results hold for ar- 
bitrary strength. The effect of strong quenched disorder 
on the force dependence of the width of transverse fluc- 
tuations cannot be simply reduced to an effective change 
in the temperature or the bending stiffness of the pure 
chain. 

We also studied the influence of uncorrelated quenched 
disorder in the transverse force along the contour of a 
WLC without spontaneous curvature under tension. In 
this case, the force-elongation curve is identical to that 
of a pure WLC with sinusoidal spontaneous curvature. 
We also find that the effect of the random force increases 
with the length of the chain. A strong random force may 
qualitatively alter the force dependence of the transverse 
fluctuations over a range of pulling forces. 

This work was supported by EPSRC via the TCM Pro- 
gramme Grant. P.B. acknowledges discussions with M. 
Hennes and K. Kroy. 
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